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Understanding the nature of superconductivity in iron-based compounds is 
essential in the development of new strategies to increase T c . Using a charge 
density analysis based on synchrotron radiation X-ray powder diffraction data, 
we found that the charge carriers only accumulated in the iron layer of the 
superconducting phase of LaFeAsOi _ X F X at low temperatures. Analysis of the 
electrostatic potential distribution revealed the concerted enhancement of the 
electronic polarization of the As ions and the carrier redistribution. This 
suggests that the enhanced electronic polarization of the As ion plays an 
important role in inducing high T c superconductivity, and that the polaron 
concept, which has been previously regarded as an untenable mechanism, 
should be reconsidered for the description of the iron-arsenide superconducting 
phase. 



1. Introduction 
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The recent discovery of new superconducting phases in iron- 
based materials (Kamihara et al, 2008; Nomura et al, 2008; 
Iimura et al , 2012) has compelled the reconsideration of the 
roles that charge, spin and lattice degrees of freedom play in 
superconductivity. For most superconducting phases observed 
at low carrier densities, it was previously thought that spin 
fluctuation plays an important role in the onset of super- 
conductivity, in which nesting between hole and electron 
pockets is the origin of Cooper-pair formation; this is 
enhanced when the FeAs tetrahedral structure forms an ideal 
tetrahedron and the pnictide height is about 1.39 A (Lee et al, 
2008; Kuroki et al, 2009). Recently, a new superconducting 
phase in the high-carrier-density region was discovered, and 
an orbital fluctuation mechanism involving the degeneration 
of Fe 3d bands was proposed as the cause of the pairing 
(Iimura et al, 2012). Thus the debate on the mechanism of 
superconductivity in iron-based materials has been reignited 
because of these various proposed pairing mechanisms. Here, 
by means of a precise charge-density study via synchrotron 
radiation X-ray powder diffraction, we report a new feature, 
electronic polarization, associated with carrier accumulation 
that is only observed in the superconducting phase of 
LaFeAsO! _ X F X . We propose that electronic polarization is 
another key factor in the onset and/or enhancement of 
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superconductivity in iron-based materials, in addition to spin 
and orbital fluctuations. 

2. Experiment 

Powder samples of LaFeAsO! _ X F X (x = 0-0.200) were 
synthesized by solid-state reaction using Fe 2 0 3 , Fe, LaAs and 
FeF 2 powders as starting materials in evacuated silica tubes. 
The fluorine contents, x, of the samples were determined by 
secondary ion-microprobe mass spectrometry. The super- 
conducting transition temperatures (T c ) were characterized as 
the onset temperatures (r c onset ) of the Meissner signals by 
superconducting quantum interference device (SQUID) 
measurement. The X-ray powder diffraction measurements on 
the powder samples of LaFeAsOx _ X F X (x = 0-0.200) were 
carried out at Powder diffraction beamline BL02B2 and 
RIKEN Materials Science beamline BL44B2 of SPring-8 to 
obtain the patterns of high statistics (~ 3 x 10 5 counts) with 
high angular resolution (d < 0.5 A). The powder samples were 
sealed in a 0.3 mm glass capillary. The temperature was 
controlled by an N 2 gas flowing system at 90-300 K and a 
helijet system (Agilent Technologies) below 90 K. The wave- 
length of incident X-rays was 0.4999 (1) A. The diffraction 
patterns were collected for 3-5 min on an imaging plate 
installed on a large Debye-Scherrer camera. Using the 
structural model obtained by Rietveld refinement, we have 
determined the precise structure parameters such as lattice 
constant, bond length, bond angle and interlayer distance 
upon F-doping. For analysis of electron density structure by 
means of Maximum Entropy Method (MEM), a total of 352 
observed structure factors were extracted after Rietveld 
refinement; the MEM analysis of each composition was 
carried out with the program ENIGMA using 
128 x 128 x 128 pixels (Sakata et al, 1990; Takata & Sakata, 
1996). The charge density obtained was converted to the 
electrostatic potential U(r) (Tanaka et al., 2006). U(r) is 
divided into the two components, that is, the nucleus charge 
component [t/ n ( r )] an d the electron charge component 
[t/ e (r)]. U n (r) is estimated by the ordinary Ewald's method, 
while t/ e (r) is directly calculated from the MEM charge 
density. The electrostatic potential obtained was visualized 
using the OpenDX program provided by IBM Visualization 
Data Explorer. 

3. Results and discussion 

The crystal structure of LaFeAsO is composed of an alter- 
nately stacked structure of FeAs 4 conducting layers and La 4 0 
blocking layers. The undoped parent compound has a struc- 
tural phase transition from a tetragonal phase to an ortho- 
rhombic phase near 150 K upon cooling, and it undergoes 
antiferromagnetic spin-density-wave (SDW) ordering (Kami- 
hara et al., 2008; Nomura et al., 2008). With F-doping on the 
oxygen sites, this structural phase transition and anti- 
ferromagnetic formation are suppressed and super- 
conductivity is exhibited (Fig. SI). To determine the structural 
dissimilarity between the superconductor and the non-super- 



conductor, we have analyzed the powder diffraction data for 
LaFeAsO! _ X F X using the MEM/Rietveld method (Sakata et 
al., 1990; Takata & Sakata, 1996). The results show that the 
fundamental structural parameters simply follow Vegard's 
rule (Denton & Ashcroft, 1991). The bond angle and the 
pnictide height of the FeAs tetrahedron, which have been 
considered as key parameters in iron-based superconducting 
materials, linearly decreased and increased, respectively, with 
F-doping. Other structural parameters such as the bond length 
and the bond angle of the LaO layer as well as the distance 
between FeAs and the LaO layer also displayed linear rela- 
tionships with F-doping (Figs. S3 and S4). The fundamental 
structural parameters were in very good agreement with 
previously reported values (Kamihara et al., 2008; Nomura et 
al. , 2008) and exhibited clear deformation of the FeAs tetra- 
hedron as a result of F-doping. There existed no clear corre- 
lation, however, between the structural parameters at the 
atomic position level and T c of this system. 

In contrast to conventional analyses of the structural 
parameters, the MEM charge density, which was obtained 
directly using the integrated intensities of the X-ray diffraction 
pattern, showed significant changes in the superconducting 
phase. Fig. 1 shows the total number and changes in the 
number of electrons at 19.5 K relative to that at 300 K for the 
Fe-centric layer and the intervening layer. Whereas the total 
number is assumed to increase monotonically as the F- 
concentration (x) increases (Fig. lb) and no temperature 
dependence is expected, the distribution of electrons in each 




(value of the n th lines = 0.5(l+«(n-l)) 



Figure 1 

F-concentration (x) dependence of the total number of electrons (b), and 
changes in the electron number within the Fe-centric layer (c) and 
intervening layer (d) at 19.5 K. (a) shows the contour map of MEM 
charge density at (100) (black contour lines) and (010) plane (grey 
contour lines) of LaFeAsO 0 95 F 0 05 . To evaluate the degree of electron 
accumulation within the Fe layer, the computing region was divided into 
two layers on the basis of the As site. An (= n — n 300 K , where n is the 
electron number within a layer) on each layer was normalized by the 
electron number (n 300 K ) at room temperature. 
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layer in the unit cell depends strongly on x and T, being 
enhanced at the values of x at which superconductivity occurs; 
at 19.5 K, the change in the number of electrons in the Fe- 
centric layer increased with x and showed a maximum around 
x = 0.100 (Fig. lc). Trends in the intervening layer were 
opposite to those in the Fe-centric layer (Fig. Id). At low 
temperatures, this fascinating feature is supposedly caused by 
the accumulation of electrons in the Fe layer and was 
confirmed by the electrostatic potential analysis of the MEM 
results (Tanaka et al., 2006; Kim et al, 2011; Fujiwara et al, 
2012). Figs. 2(a)-(c) show two-dimensional maps with equi- 
potential lines. Unlike the MEM charge density, which may be 
affected by hybridization, the electrostatic potential allows the 
volume (V EP ) occupied by each atom to be estimated (Fuji- 
wara et al, 2012). The three-dimensional volume is formed by 
traces of the local minimum point of the electrostatic poten- 
tial. According to the electrostatic potential formula, the 
boundary of V EP is placed in a quasi-neutral point; an increase 
in the number of electrons causes a shrinkage of V EP to 
maintain neutral character in the boundary of electrostatic 
potential. We focused on variations in the V EP of the Fe atom 
(V E p), as the behavior of the Fe atom is thought to be closely 
linked to superconductivity in this system. At 19.5 K, V E p 
shrank by over 7% at x = 0.050, but it then increased at x = 
0.150 (Fig. 2e). This large shrinkage cannot be attributed to a 
decrease in the lattice volume alone, which resulted in at most 
a 1.3% decrease at x = 0.200 (Table SI). The variation in V E p, 
therefore, reflects that electrons are effectively concentrated 
in the Fe layer in the 0.050 < x < 0.150 region at this low 
temperature. 

Along with this electron charge redistribution in the Fe 
layer, some notable features have been observed in the 
variations of the local electrostatic potential. At x = 0.125, the 
equipotential lines showed distinct differences in the upper 
and lower parts of the As ion, compared with isotropic 
variations at x = 0 and 0.200; the upper equipotential lines 
show a gradual slope with smooth curves, in comparison to the 
closely spaced lower equipotential lines. These effects lead to 
an extremely asymmetrical distribution of the electrostatic 
potential around the As ion, which arises from a separated 
electrostatic potential distribution on the As ion (Fig. S5). 
Such a separated electrostatic potential distribution on an 
atom was observed from electrostatic potential imaging of 
PbTi0 3 by Tanaka et al. (2006), who reported that the specific 
separation of the electrostatic potential distribution is a 
signature of polarization. Enhancement of the electronic 
polarization induces a large change in the electrostatic 
potential distribution because of its influence on the electro- 
static interactions between electrons distributed in that space 
and those on the As ion. In order to estimate the degree to 
which the electronic polarization of the As ion was enhanced, 
we evaluated the asymmetric ratio, Pi/P 2 , where P 1 (P 2 ) is the 
minimum value of electrostatic potential between the As ion 
and lower (upper) La ion. In the superconducting phase (x = 
0.125), P 1 and P 2 were respectively -25.24 and -20.06 V 
(P\IP 2 = 1.258) (Fig. 2d). A noticeable difference was not 
observed at x = 0 (P 1 /P 2 = 1.002) and x = 0.200 (P 1 /P 2 = 1.096), 



which were not the superconducting phase (Fig. 2f). It is worth 
noting that the electronic polarization of the As ion is signif- 
icantly enhanced in the superconducting phase, which was 
clearly shown by the asymmetric contour curves and the 
resultant asymmetrical spatial distribution of the electrostatic 
potential around the As ion. 

Both y E p and the asymmetric ratio of the electrostatic 
potential on the As ion vary not only with F-concentration (x) 
but also with temperature, as shown in Fig. 3. At room 
temperature for x = 0.125, V E p was reduced significantly by 
13%, indicating that surplus electrons from F-doping are 
confined in the Fe layer at 300 K (Fig. 3a). However, no 
further shrinkage of V E p was observed at x = 0.200, which 
suggests that the number of electrons that can accumulate in 
the Fe layer reaches its maximum for x = 0.125 at 300 K, 
although the total number of electrons continues to increase 
as x increases. After cooling to 19.5 K, however, V E p for x = 
0.125 shrank by a further 7%, indicating further electron 
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Figure 2 

Electrostatic potential map with equipotential lines (a)-(c), variation of 
the electrostatic potential volume (V^p) of the Fe atom (e) and 
asymmetric ratio (Pi/P 2 ) of electrostatic potential distribution around 
the As ion (/) as a function of F-content (x) at 19.5 K. T c are plotted in 
(g). The dashed curve around Fe in (a)-(c) is the trace of the local 
minimum points of electrostatic potential around the Fe atom. The 
minimum points, P 1 (P2), of electrostatic potential distribution around the 
As ion are shown in (d) which is the one-dimensional plot of the dashed 
line crossing the center position of the La and the As ion. T° nset and T c s in 
(g) are defined by the onset temperatures of the Meissner signals and the 
temperatures of the structural phase transition from a tetragonal phase to 
an orthorhombic phase. 
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accumulation. Although V^> for x = 0 also shrank with 
decreasing temperature, at 19.5 K it had not reached the value 
for x = 0.125 at 300 K; an insufficient number of electrons 
accumulated for superconductivity. Vep at jc = 0.200, on the 
other hand, increased by ~ 7%, showing that electron 
numbers in the Fe layer decreased with cooling. To account for 
the large contraction of Vep at jc = 0.125, an additional effect, 
which can outweigh the less-than-1% shrinkage of the lattice 
volume, is required at low temperature. As well as the marked 
contraction of V^>, P\IPi for the As ion increased consider- 
ably at x = 0.125, for which superconductivity was observed 
under 26 K (Fig. 3b). As the temperature decreased, P 1 /P 2 
scarcely changed for x = 0. P\IP 2 of both doped samples, on the 
other hand, began to increase near 90 K, and P\\P 2 for x = 
0.125 rose significantly to a value of 1.26 until the temperature 
had fallen to 19.5 K, whereas that for x = 0.200 increased no 
further than 1.1 at 42 K. 

Based on V^> and P\IP 2 , we can understand the effect of 
doping on the emergence of superconductivity as follows. At 
300 K surplus electrons from F-doping are confined to the Fe 
layer and the electronic polarization of the As ion is not 
changed. As temperature decreases from near 70 K, the 
electronic polarization of the As ion is only significantly 
enhanced in the superconducting phase, and the accumulation 
of electrons in the Fe layer is accelerated. Consequently, our 
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Figure 3 

Temperature dependence of the electrostatic potential volume (V^f P ) of 
the Fe atom (a) and the degree of asymmetric space distribution (P 1 /P 2 ) 
around the As ion (b) at x = 0.125 (red circles) and x = 0.200 (blue circles). 
Black circles are the data for x = 0 at 19.5 and 300 K. 



results suggest that the electronic polarization of the As ion is 
one of the crucial factors in inducing superconductivity in 
LaFeAsOi _ X F X . We postulate that this electronic polarization 
is the average result, as observed by the steady-state X-ray 
diffraction, of polaronic behavior, which is a dynamic motion. 
Until now, electron-phonon coupling has been regarded as an 
untenable mechnism in high T c iron-based superconductivity, 
because a T c value beyond the McMillan limitation (McMillan, 
1968) could not be explained by such a pairing mechanism 
(Kamihara et al, 2008; Nomura et al, 2008; Iimura et al, 2012). 
In spite of this weakness, an additional effect related to 
polarization has been consistently proposed in iron super- 
conductivity. Recently Sawatzky et al. has suggested a 
bipolaronic concept (Sawatzky et al, 2009; Berciu et al, 2009), 
in which the fluctuation of electrons on the Fe ion leads to the 
electronic polarization of the As ion; the polarized electron 
clouds interact attractively, and thereby propagate effective 
charge carriers into the Fe layer. According to the theory, the 
As 4p-orbital is polarized toward the Fe atom by residual 
electrons on the Fe atom. As shown in Figs. S5 and S6, the 
electronic polarization of the As ion is, as an average of the 
dynamic fluctuation, in the direction of the c-axis. This is very 
likely indicative of a polaron or bipolaron; the in-plane 
component of the polaron was canceled out in the time- 
averaged XRD results, and the polarization induced by the 
bipolaron was observed along the odirection. Although only 
three cases have been computed and the many-body problem 
still remains, the direction of electronic polarization being 
aligned toward the c-axis, which has been visualized by elec- 
trostatic potential imaging, is reasonable evidence for the 
existence of polaron/bipolaron, in addition to the enhanced 
electronic polarization only occurring in the superconducting 
phase. As well as the specific direction and emergence of the 
electronic polarization of the As ion in the superconducting 
phase, the onset temperature of electronic polarization also 
suggests the existence of polaron. Recently, Zhang et al (2010) 
measured the Fe— As mean-square bond displacements 
((TAs-Fe) m LaFeAsOi _ X F X (x = 0 and 0.07) using X-ray 
absorption fine-structure techniques. Generally, the 0As-Fe 
values decreased as the temperature decreased due to the 
suppressed thermal vibration of the atom. However, Oas— Fe 
for x = 0.07, for which superconductivity appears at 26 K, 
decreased with cooling at temperatures above 70 K, but it then 
increased abruptly near 70 K. This anomaly reached a 
maximum near 30 K, while cr As _ Fe for x = 0 decreased 
continuously as the temperature decreased, and the same 
value remained almost constant at temperatures below 150 K. 
It was proposed that the abnormal enhancement of cr As _ Fe was 
due to polaron formation in the superconducting phase of 
LaFeAsO0.93F0.07- The polarization of the As ion seen here, 
shown as P\IP 2 in Fig. 3(b), is in good agreement with this 
result. That is, the asymmetric spatial distribution can be 
related to the distorted electron clouds in the FeAs layers, and 
implies that superconductivity in this system can be ascribed 
to polarons caused by electronic polarization of the As ion. 
The good correspondence between the carrier accumulation in 
the Fe layer, the electronic polarization of the As ion, and the 
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appearance of superconductivity suggests that the electronic 
polarization of the As ion and carrier accumulation are 
cooperative factors in the emergence of iron-arsenide super- 
conductivity. 



4. Summary 

In summary, we have visualized the electronic polarization of 
the As ion within LaFeAsOi _ X F X via MEM charge density 
and electrostatic potential methods using synchrotron radia- 
tion X-ray powder diffraction data. It was found that 
enhanced electronic polarization of the As ion is one of the 
most important factors, along with the accumulation of 
carriers in the Fe layer, for the onset of high T c super- 
conductivity. Moreover, we have demonstrated that the 
polaron mechanism mediated by electronic polarization of the 
As ion together with characteristics of the Fe ion should not be 
overlooked in the analysis of iron-arsenide superconductivity. 
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